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ABSTRACT 

We present the first interferometric 12 CO(J = 3-2) observations (beam size 
of 3'.'9 x l'/6 or 160 pc x 65 pc) with the Submillimeter Array (SMA) toward 
the center of the Seyfert 2 galaxy M51. The image shows a strong concentration 
at the nucleus and weak emission from the spiral arm to the northwest. The 
integrated intensity of the central component in 12 CO(J = 3 - 2) is almost twice 
as high as that in 12 CO(J = 1-0), indicating that the molecular gas within 
an ~ 80 pc radius of the nucleus is warm (> 100 K) and dense (~ 10 4 cm -3 ). 
Similar intensity ratios are seen in shocked regions in our Galaxy, suggesting that 
these gas properties may be related to AGN or starburst activity. The central 
component shows a linear velocity gradient (~ 1.4 km s _1 pc -1 ) perpendicular 
to the radio continuum jet, similar to that seen in previous observations and 
interpreted as a circumnuclear molecular disk/torus around the Seyfert 2 nucleus. 
In addition, we identify a linear velocity gradient (0.7 km s _1 pc -1 ) along the jet. 
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Judging from the energetics, the velocity gradient can be explained by supernova 
explosions or energy and momentum transfer from the jet to the molecular gas 
via interaction, which is consistent with the high intensity ratio. 

Subject headings: galaxies: individual (M51, NGC 5194), galaxies: ISM, galaxies: 
nuclei, galaxies: Seyfert 

1. INTRODUCTION 

According to the unified scheme as originally suggested by Antonucci & Miller (1985), 
two types (Types 1 & 2) of Active Galactic Nuclei (AGNs) originate from the same object (a 
supermassive black hole surrounded by an optically-thick torus) that is viewed at different 
angles. This hypothesis has stimulated intensive millimeter-wave interferometric observa- 
tions to search for molecular tori toward AGNs. For example, images of the prototypical 
Seyfert 2 galaxy NGC 1068 show a ~ 100 pc disk or torus of molecular gas around its AGN, 
with a rotational axis parallel to the axis of its radio continuum jet (Jackson et al. 1993; 
Tacconi et al. 1994; Heifer & Blitz 1995). Sub-arcsecond 12 CO(J = 2-1) observations sug- 
gest that this torus/disk is warped and viewed nearly edge-on at a radius of less than 110 pc 
(Schinnerer et al. 2000). On the other hand, recent surveys of Seyfert galaxies suggest that 
the circumnuclear molecular gas has a wide variety of distributions, and does not always ex- 
hibit a centrally-peaked component (e.g., Kohno et al. 2001). It may be that only a fraction 
of AGNs have molecular tori detectable in current observations, making detailed studies of 
such examples especially important. 

M51 (NGC 5194), a grand-design spiral galaxy at a distance D ~ 8.4 Mpc (Feldmeier, 
Ciardullo, & Jacoby 1997), has a Type 2 AGN (Ho, Filippenko, & Sargent 1997). As in 
the case of NGC 1068, Kohno et al. (1996) showed that the AGN of M51 is surrounded by 
a ~ 100 pc disk/torus of dense molecular gas with a rotational axis parallel to the axis of 
the central radio jet (Ford et al. 1985). The estimated column density of the molecular gas 
exceeds 10 23 cm' 2 , which is comparable with the amount of obscuring material inferred from 
X-ray observations (e.g., Terashima & Wilson 2001). 13 CO and HCN J = 1 - observations 
suggest that the circumnuclear molecular gas is dense (~ 10 5 cm~ 3 ) and warm (> 300 K; 
Matsushita et al. 1998, 1999). These conditions suggest that the nuclear disk/torus can be 
better studied in higher-J transitions, such as 12 CO J = 3 - 2 (the critical density for this 
transition is ~ 5 x 10 4 cm -3 , and the energy of the J = 3 level corresponds to 33 K ). 

Here, we report observations of the central region of M51 with the Submillimeter Array 7 
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(SMA; Ho, Moran, & Lo 2004) in the 12 CO(J = 3-2) line. As anticipated, the central 
molecular torus/disk emits stronger in the J = 3 - 2 than the J = 1 - 0, but we also detected 
emission from a spiral arm. We point out that this torus/disk exhibits a velocity gradient 
not only along its major axis, but also along its minor axis parallel to the radio jet. 

2. OBSERVATIONS 

The submillimeter-wave aperture synthesis images toward the center of M51 were ob- 
tained in the 12 CO(J = 3-2) line (rest frequency = 345.796 GHz) with the SMA. The 
observations were made on 2003 June 12 and 14 during the testing and commissioning phase 
when a total of four 6 m antennas were available. The SMA correlator had a 0.96 GHz band- 
width at the time, and was configured to have a frequency resolution of 0.8125 MHz. We 
observed Mars and Uranus, respectively, for bandpass and flux calibration, and J1310+323 
every 20 minutes for amplitude and phase calibration. The uncertainty in the absolute flux 
scale is estimated to be ~ 20%. 

We calibrated the data using the OVRO software package MIR, which was modified 
for the SMA. The images were CLEANed using the NRAO software package AIPS, and the 
resultant angular resolution was 3'.'9 x 1'.'6 (160 pc x 65 pc) at a position angle (P.A.) of 146° 
with natural weighting. We made channel maps with 42 channel binning, which corresponds 
to a velocity resolution of 29.6 km s _1 . The typical rms noise level of the channel maps 
was 70.4 mJy beam™ 1 , which corresponds to T sySj DSB ~ 700 K. The half-power width of the 
primary beam at 345 GHz is 36" (1.5 kpc). We did not detect any continuum emission at 
the rms noise level attained of 18 mJy beam™ 1 . 

3. RESULTS 

Channel and integrated intensity maps in 12 CO (J = 3 - 2) as shown in Figures 1 and 
2 show a strong peak at the nucleus at redshifted velocities (with respect to the systemic 
velocity of 472 km s _1 ; Scoville & Young 1983; Tully 1974). The highest velocity of the 
compact nuclear emission is ~ 616 km s™ 1 (as can be seen also in the position-velocity 
diagram of Fig. 3), comparable with the highest velocity seen in the redshifted wing of the 
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12 CO(J = 3-2) spectrum taken with the JCMT (Matsushita et al. 1999). Our maps also 
show weak emission from the spiral arm to the northwest of the nucleus (hereafter, the 
northwestern arm) at blueshifted velocities. The spiral arm to the south-east (hereafter, the 
southeastern arm) was not detected. The previously published single-dish 12 CO(J = 3-2) 
maps (Wielebinski et al. 1999) also show the intensity asymmetry, which is stronger emission 
from the northwestern than from the southeastern arm. 

The compact nuclear emission also appears at redshifted velocities in interferometric 
12 CO(J = 1-0) observations by Aalto et al. (1999) and Sakamoto et al. (1999). Similarly, 
interferometric 12 CO(J = 2-1) observations by Scoville et al. (1998) also show predominantly 
redshifted emission, concentrated 1" west of the nucleus, but in addition relatively weak 
blueshifted emission east of the nucleus (see their Figs. 1 and 2). The integrated intensity of 
their blueshifted emission is only 18% of that of the redshifted emission, with a deconvolved 
size for the blueshifted emission of less than 0"5. In our map, the redshifted emission was 
detected at a peak integrated intensity of 15a at a beam size of 3" 9 x 1"6. Assuming the 
redshifted and blueshifted emission in J = 3 - 2 has the same intensity ratio as that in J = 
2-1, the blueshifted emission is expected to have a signal-to-noise ratio of < 2a. We rule 
out the possibility that the systemic velocity of the circumnuclear molecular gas is different 
from that of the host galaxy (472 km s _1 ). Interferometric HCN(J = 1-0) observations by 
Kohno et al. (1996) show a compact nuclear component with comparable bright redshifted 
and blueshifted emission symmetrically placed either side of the nucleus, interpreted as an 
approximately edge-on rotating disk or torus. 

The missing flux in our map was estimated by comparing our data with previous single- 
dish observations. The flux within the central 13" and 24" region in our 12 CO(J = 3-2) 
map (corrected for primary beam attenuation) is 13.8 K km s _1 (228 Jy beam -1 km s _1 ) 
and 8.8 K km s _1 (499 Jy beam -1 km s _1 ) respectively. The single-dish flux are 43.1±2.0 
K km s" 1 within 13" (Matsushita et al. 1999) and 47±7 K km s" 1 within 24" (Dumke et 
al. 2001). We therefore recovered only 32% and 19% of the single-dish flux for the central 
13" and 24" respectively, suggesting that there exists extended 12 CO(J = 3-2) emission in 
the central region. Indeed, the single-dish 12 CO(J = 3-2) maps show very similar image 
with the single-dish J = 2- lorl-0 maps, namely, even J = 3 - 2 emission is rich in 
extended structures (Wielebinski et al. 1999; Dumke et al. 2001), which support our result. 
The missing flux of ~ 70 — 80% is, however, larger than that in the past extragalactic 12 CO(J 
= 1-0) observations of < 50% (e.g., Sakamoto et al. 1999). The possibility of calibration 
errors in single-dish data should be low, because the above values are consistent with each 
other, and also consistent with the independent HHT (49 ±2 K km s -1 ; Mauersberger et 
al. 1999) and the CSO (48.3 K km s" 1 with 24" beam; Bash et al. 1990) observations. An 
underestimate of the absolute flux in our data is possible, but this would suggest more 
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extreme conditions at the central region of M51 (see the next paragraph and section). 

We compared the intensities of the nuclear molecular gas components using our 12 CO(J 
= 3-2) data with the interferometric J = 1 - data (Sakamoto et al. 1999). First, we 
convolved our J = 3 - 2 image to the same resolution as that in J = 1 - (4'.'2 x 3'.'4). The 
J = 3 - 2 and 1-0 line intensities at the nucleus averaged over this beam size (170 pc x 
140 pc) are 91.7 ± 4.8 K km s" 1 and 49.2 ± 5.2 K km s -1 , respectively (temperatures are 
in the Rayleigh- Jeans approximation, and the uncertainties are ±1<t). The J = 3 - 2 line is 
therefore stronger than the J = 1 - line, with a 12 CO(J = 3 - 2)/(J = 1-0) intensity ratio, 
i?3i, of 1.9 ± 0.2. This line ratio may be a lower limit, as the J = 1 - data have shorter 
baselines than the J = 3 - 2 data. To remedy this, we truncated the inner Mf-plane of the 
J = 1 - data to match that in the J = 3 - 2 data, and recomputed the line intensities 
following the above method. The resultant R S i at the nucleus averaged over this beam size 
(3'.'9 x 2'.'6 or 160 pc x 110 pc) is ~ 5.2 ± 1.7, where the larger uncertainty reflects the noisier 
J = 1 - image due to the inner uv truncation. In the following discussion, we adopt R31 
of 1.9 ± 0.2, keeping in mind that this is probably a lower limit. 

We made position-velocity (PV) diagrams along and perpendicular to the axis of the 
radio continuum jet (P. A. of ~ 165°) from the nucleus (Fig. 3). To compute the magnitude 
of the velocity gradient, we first searched for the strongest peak (exceeding 2a) within a 5" 
radius of the nucleus at each velocity bin (width of 14.8 km s" 1 ) in the PV diagram. We 
then applied a linear fit to these peaks weighted by their respective intensities. In this way, 
we inferred a linear velocity gradient perpendicular to the jet axis of 1.39 ±0.01 km s^ 1 pc -1 , 
and parallel to the jet axis of 0.65 ± 0.01 km s _1 pc -1 . The sense of the velocity gradient 
perpendicular to the jet (Fig. 3a) is the same as that found in the past interferometric 
observations (Kohno et al. 1996; Scoville et al. 1998); the magnitude of the velocity gradient 
of 1.39 km s _1 pc -1 is comparable with that estimated from the interferometric HCN(J = 
1-0) observations of 1 — 2 km s _1 pc -1 (Matsushita et al. 1998). This velocity gradient 
is attributed to a nearly edge-on rotating molecular disk or torus around the nucleus. The 
velocity gradient along the jet has not been previously discussed, but can also be seen in the 
channel/velocity maps in both CO J = 1 - and 2-1 (see Fig. 2 of Scoville et al. 1998 
and Fig. lq of Sakamoto et al. 1999). Based on these arguments, we believe the existence 
of a velocity gradient parallel in addition to that perpendicular to the radio jet to be real, 
although we note that the circumnuclear emission is poorly resolved and hence the actual 
magnitude of the velocity gradients may need to be revised when observations at higher 
angular resolutions become available. 
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4. DISCUSSION 

The i?3i within an ~ 80 pc radius of the nucleus of M51 is 1.9 ± 0.2, more than twice as 
high as the ratios measured by single-dish telescopes at lower spatial resolutions; line ratios 
observed with the beam sizes of ~ 14" — 24" (~ 600 — 1000 pc) are 0.5 — 0.8 (Matsushita 
et al. 1999; Mauersberger et al. 1999; Wielebinski et al. 1999). The physical properties of 
the circumnuclear molecular gas must therefore be considerably different from that in the 
disk of the galaxy. We estimated the physical conditions of the circumnuclear molecular gas 
using the Large- Velocity-Gradient (LVG) approximation (Goldreich & Kwan 1974; Scoville 
& Solomon 1974) assuming a one-zone model. The collision rates for CO in the temperature 
range 10 — 250 K were taken from Flower & Launay (1985) and 500 — 2000 K from McKee et 
al. (1982). The velocity gradient is estimated to be about 1.39 km s~ 4 pc -1 (Sect. 3). This 
value is consistent with the gradient caused by the internal turbulence (~ 2 km s -1 pc -1 ; 
for a more detailed discussion, see Matsushita et al. 1998). Using the 'standard' relative 
abundance [ 12 CO]/[H 2 ] = 5 x 10~ 5 , we find the H 2 number density n H2 > 10 4 cm -3 and 
the kinetic temperature T k > 500 ± 150 K, with the uncertainties reflecting the uncertainty 
in i? 3 i. Even for an order of magnitude lower relative abundance or an order of magnitude 
higher velocity gradient, this model gives nn 2 > 10 4 cm~ 3 and Tk > 95 ± 15 K. We therefore 
conclude that the circumnuclear molecular gas is denser and warmer than the molecular gas 
in the galactic disk. This is consistent with the results of the interfero metric 13 CO and HCN 
J = 1 - observations, which also suggest dense (riH 2 ~ 10 5 cm -3 ) and warm (T^ > 300 K) 
molecular gas based on LVG calculations (Matsushita et al. 1998, 1999). 

There are several possible reasons for the rather extreme conditions of the molecular 
gas. One is interaction with the radio jet. Molecular outflows from star forming regions in 
our Galaxy show R 31 greater than unity (e.g., Richardson et al. 1985; Hirano & Taniguchi 
2001), similar to that in the center of M51. The strongest 12 CO(J = 3-2) emission in our 
channel maps comes from 1" north-west of the nucleus (484 km s -1 in Fig 1), and this also 
is a region with high [N II] /Ha intensity ratio (~ 4). The high [N II] /Ha ratio suggests that 
relativistic particle heating is dominant (Cecil 1988), which suggests an interaction between 
the jet and molecular gas. Recent high resolution H2O maser observations reveal that the 
masers are located close to the nucleus, and are possibly associated with the jet (Hagiwara 
et al. 2001). The velocity range of the H2O masers of 538 — 592 km s _1 is similar to that 
of the redshifted nuclear emission in our data (see Fig. 1 and 3), which also supports the 
idea that the molecular gas is interacting with or entrained by the jet. Shock/compression 
by supernova (SN) explosions is another possibility, since shocked molecular gas around SN 
remnants in our Galaxy also shows high intensity ratios (e.g., Arikawa et al. 1999). 

The linear velocity gradient along the jet axis is an interesting phenomenon, and it is 
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worth considering what may cause this gradient. We first calculate the mass, momentum, 
and kinetic energy of the molecular gas within the velocity range of 515 - 615 km s _1 , which 
is away from the systemic velocity. The CO-to-H 2 conversion factor can be estimated using 
the LVG approximation mentioned above, and the result is ~ 1 x 10 20 cm -2 (K km s -1 ) -1 . 
The mass is calculated using this conversion factor and the intensities of molecular gas. In 
each channel map within the velocity range, the velocity of the molecular gas is defined as 
the velocity difference from the systemic velocity, and we assume that the molecular gas 
is moving along the line of sight. Thus the estimated values are lower limits. Adding all 
channel map information, we estimated the mass, momentum, and kinetic energy as 1 x 10 6 
M , 2 x 10 46 g cm s" 1 , and 8 x 10 52 erg, respectively. 

Using these values, we discuss the possible causes of the velocity gradient along the jet 
axis. One is that the molecular gas away from the systemic velocity is entrained by the jet, 
namely, relativistic particles. The minimum energy of the jet is 6.9 x 10 51 erg (Crane & van 
der Hulst 1992), and the momentum is 2 x 10 41 g cm s _1 , if we assume that the velocity of 
the relativistic particles is 90% of the speed of light. Crane & van der Hulst (1992) estimated 
the energy of the jet assuming the total energy of the jet plasma is 100 times larger than 
that of electrons. This number can be larger, but it should be smaller than the mass ratio of 
proton and electron of about 1800. The actual geometry of the jet and the molecular gas is 
unclear, but should be similar if molecular gas is entrained by the jet. Therefore the energy 
of the jet may be similar to that of the molecular gas, but the momentum is several orders 
of magnitude different. This problem is also seen in the study of young stellar objects, 
and still under discussion (e.g., Richer et al. 2000). Another possibility is that since the 
momentum and energy can be transferred from one to another by interaction. If those of the 
molecular gas along the jet are continuously input by the jet via interaction, the momentum 
and energy differences between the molecular gas and the jet can be explained. Indeed, the 
high [N II] /Ha ratio (an evidence of relativistic heating) and the H 2 masers (an evidence 
of shocked molecular gas) are observed, as mentioned above. These observations support the 
idea that the energy and momentum of the molecular gas along the jet are transferred from 
the radio jet via interaction. 

There are several possible explanations to create the velocity gradient along the jet 
without considering interaction with the radio jet. Supernova (SN) explosions can be one of 
the possible causes, based on the similarity of R 3 i as mentioned above. However, 10 2 — 10 4 
SNe are needed to explain the energy of the molecular gas, if we assume that one SN releases 
10 51 erg and that the energy transfer efficiency is at most 20% (McCray & Kafatos 1987). 
Alternatively, the velocity gradient along the jet may be due to a warped disk around the 
AGN (e.g., Schinnerer et al. 2000). In this case, not all the gas is on a plane of a molecular 
gas disk. Thus, even if most of the gas show velocity gradient perpendicular to the jet, some 
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of the gas show velocity gradient along the jet. Future SMA observations with higher spatial 
resolution will provide more detailed information about the molecular gas around the AGN. 

We thank P.T.P. Ho, J.M. Moran, and M.J. Cai for valuable comments. We also thank 
all the past and present SMA staffs for designing, constructing, and supporting the SMA. 
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Fig. 1. — Channel maps of the 12 CO(J = 3-2) line of the central region of M51. The synthesized 
beam (3'.'9 x l'/6 or 160 pc x 65 pc) is shown at the lower-left corner of the first channel map. LSR 
velocities are shown at the bottom-right corner of each map. The systemic velocity is 472 km s _1 
(Scoville & Young 1983; Tully 1974). The cross and the 36" (1.5 kpc) diameter circle in each map 
are the galactic nucleus determined from the 6 cm radio continuum peak (Ford et al. 1985) and the 
half-power width of the primary beam. The contour levels are —6, —3, 3, 6, 9, 12, and 15a, where 
la = 70.4 mJy beam -1 (= 0.115 K). Primary beam correction is not applied. 
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Fig. 2. — Integrated intensity map of the 12 CO(J = 3-2) line. The gray scale and the 
synthesized beam are shown at the top and the bottom-left corner, respectively. The cross 
and the circle are the same as in Fig 1. The contour levels are 2, 3, 4, 6, 9, 12 and 15a, where 
la = 5.89 Jy beam -1 km s _1 (= 9.63 K km s _1 ). Primary beam correction is not applied. 
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Fig. 3. — Position-velocity diagrams (a) along and (b) perpendicular to the radio jet (P.A. 
~ 165°). The velocity resolution is 14.8 km s" 1 . Positions are the offsets from the nucleus. 
The contour levels are —2,2,4,6,8,10 and 12a, where la = 76.5 mJy beam -1 for (a) and 
47.8 mJy beam" 1 for (b). Dotted horizontal lines are the systemic velocity. Dashed lines 
show the results of intensity-weighted linear fitting. 



